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Abstract 

We report the detection of a time variable OVII line emission in a deep 100 ks 
Suzaku X-ray Imaging Spectrometer spectrum of the Galactic Ridge X-ray emission. 
The observed line intensity is too strong (11±2 line unit or photon cm~^ s~^ str~^) to 
be emitted inside the heavily obscured Galactic disk. It showed a factor of two time 
variation which shows a significant (~4(t) correlation with the solar wind O^"*" ion flux. 
The high line intensity and the good time correlation with the solar wind strongly 
suggests that it originated from geocoronal solar wind charge exchange emission. We 
discuss the X-ray line intensity considering a line of sight direction and also theoretical 
distribution models of the neutral hydrogen and solar wind around the Earth. Our 
results indicate that X-ray observations of geocoronal solar wind charge exchange 
emission can be used to constrain these models. 

Key words: X-ray: diffuse background — Sun: solar wind — Sun: solar- 
terrestrial relations — Earth 



1. Introduction 

The solar wind charge exchange (SWCX) emission has been recognized to reproduce 
X-ray emission in many solar system environments such as comets, geocorona, heliosphere, and 
planetary atmospheres (see Bhardwaj et al. 2007 for review). When an ion in the solar wind 
interacts with a neutral atom, it strips an electron(s) from the atom and X-ray or ultra-violet 
photon(s) are released as the electron relaxes into the ground state. Therefore, the SWCX 
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spectrum is characterized by strong emission lines from highly ionized atoms such as O or 
C^'^. It can be distinguished from ordinary thermally equilibrium plasma emission by strong 
emission lines from large principle quantum number states such as n=3-6, and strong forbidden 
lines. The classical over-the barrier model (Ryufuku et al. 1980; Mann et al. 1981) can explain 
the large n line feature. Because the heavy solar wind ions are highly charged, their ionization 
potentials are larger than those of the neutral target. Therefore, when the electron moves from 
the target to the ion, it transfers into an excited state of the ion whose potential most closely 
matches that of the neutral. Consequently, the electron will enter into a high n state and then 
emit a strong emission line by the transition. The strong forbidden line feature is explained by 
a statistical-distribution model (Krasnopolsky et al. 2004). 

Among many classes of SWCX objects in the solar system, the geocorona is now consid- 
ered as an important contamination. The low density neutral hydrogen around the Earth or the 
geocorona emits X-ray emission lines such as O VII, O VIII, and C V Kq, by the SWCX. When we 
would like to observe soft X-rays from objects such as clusters of galaxies, this component can 
be a contamination source. At the same time, SWCX X-ray emission can provide a diagnostic 
of the spatial distribution of the geocorona and solar wind near the Earth. 

During the ROSAT all-sky survey, Snowden et al. (1994) found unknown long term 
enhancements (LTE) in the soft X-ray background. Motivated by the discovery of cometary 
X-rays (Lisse et al. 1996; Dennerl et al. 1997), Cox (1998) suggested that the SWCX due to 
geocorona can explain a part of the LTE. Robertson & Cravens (2003) constructed a model to 
simulate a spatial distribution of the geocoronal SWCX emission using the Earth's exosphere 
and solar wind distributions. 

Observational results, which supports these predictions, have been taken with Chandra, 
XMM-Newton, and Suzaku. Wargelin et al. (2004) found a signature of an oxygen emission 
line in the Chandra Moon observation. Because the Moon occults any background sources, it 
must have originated from the geocoronal SWCX. Snowden et al. (2004) discovered a significant 
enhancement of the soft X-ray background during one of XMM-Newton's observations of the 
Hubble Deep Field-North. The emission showed lines from O, Ne and Mg and also was time 
variable. This phenomena was nicely explained by XMM-Newton's changing viewing geometry 
with respect to the geocorona, whose SWCX is predicted to have a specific distribution around 
the Earth (Robertson & Cravens 2003). 

Carter & Sembay (2008) systematically investigated XMM-Newton archival data focus- 
ing on time variability of the soft X-ray background. Using solar wind data and considering 
the XMM-Newton's line of sight, similar to Snowden et al. (2004), they found that, in several 
cases, the SWCX enhancement occurred when XMM-Newton was on the sub-solar side of the 
magnetosheath which is predicted by Robertson & Cravens (2003) to be a high SWCX flux 
region. The SWCX-like emission was found for other cases when XMM-Newton had a line 
of sight which did not intersect the high flux region. They considered that these probably 
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arise from a non geocoronal origin such as a coronal mass ejection (CME) passing through the 
hehosphere. 

CMEs have been used to explain diffuse X-ray background in other X-ray observations. 
For example, Carter & Sembay (2010) recently reported a diffuse variable X-ray emission de- 
tected with XMM-Newton. They analyzed its spectrum and light curve in combination with 
the solar wind proton data, and concluded that this is associated with a CME which interacted 
with neutrals in the Earth's exosphere. 

Thanks to a good energy response of the CCDs onboard Suzaku, Fujimoto et al. (2007) 
discovered a firm line of evidence for the geocoronal SWCX emission in the north ecliptic pole 
observation (RA= 272.800°, Dec= 66.000°, J2000). Although the observation was originally 
conducted to take the sky X-ray background, they found a variability in X-ray CCD light 
curves, which apparently correlates with the solar wind proton flux. They divided the data 
into two periods, i.e., the stable and flare durations and successfully found enhancement of 
emission lines such as the CVI 4p to Is line and the O VII forbidden line, in the flare spectrum. 
To explain short-term variations of less than 1 hr, they introduced a new parameter, the point 
where the line of sight encounters the magnetosheath (see figure 8 in Fujimoto et al. 2007). 
They concluded that the short-term X-ray variability have ant i- correlation with the distance 
to the magnetosheath. 

In this paper, we studied a possible SWCX oxygen line emission discovered in the long 
Suzaku observation of the Galactic ridge (RA= 281.000°, Dec= -4.070°, J2000), (/ = 28.46°, b = 
—0.20°). In Ebisawa et al. (2008), we found strong O VII emission which is difficult to attribute 
to the Galactic ridge X-ray emission. Although we thought that it can be the SWCX emission, 
investigations are not enough to conclude that. In this paper, to confirm this inspection, we 
conducted a timing analysis by combining the Suzaku data with the solar wind data. While 
the north ecliptic pole observation was in the direction of the north pole, the Galactic ridge is 
near the celestial equator and the sub-solar side of the magnetosheath, where the geocoronal 
SWCX is expected to be strong. The good statistics 100 ks Suzaku data provides us with a 
good opportunity to characterize the geocoronal SWCX emission. 

2. Observation 

The Suzaku observation of the Galactic Ridge direction was conducted from 2005 
October 28, 02:24 to October 30, 21:30 for 100 ks. October 28 and 30 corresponds to the 
Day of Year (DOY) in 2005 of 301 and 303, respectively. Suzaku (Mitsuda et al. 2007) is 
the fifth Japanese X-ray astronomical satellite which carries four X-ray CCDs (X-ray Imaging 
Spectrometer, XIS: Koyama et al. 2007). Due to the low-earth orbit and the large effective 
area, the back-illuminated type CCD (XISl) has one of the lowest relative particle backgrounds 
among all X-ray CCDs in currently available X-ray observatories. Furthermore, the XISl has 
good energy resolution and superior low-energy response with negligible low-energy tails. These 
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two characteristics make Suzaku ideal for studying diffuse soft X-ray emission like the geocoro- 
nal SWCX. 

Figure 1 shows the XIS spectra in the low energy band taken from Ebisawa et al. 
(2008). Among many emission lines such as FeXVII (0.829±0.004 keV), NeX (1.021±0.002 
keV), MgXI (1.342±0.003 keV), and SiXIII (1.843±0.003 keV), we can notice a strong OVII 
line at 0.560±0.03 keV. These energies are the best fit values in table 2 of Ebisawa et al. (2008). 
The line center energy of OviI is more likely the forbidden line (561 eV) than the resonance 
line (574 eV). Its intensity is too strong to arise from the Galactic ridge X-ray emission which 
must suffer from the interstellar hydrogen column density of A^h ~ 1 x 10^^ cm^. The flux is 
(320±50)xl0~^ photons s~^ cm~^ deg~^, which corresponds to 11±2 line units (LU, photons 
s~^ cm~^ str~^). This is even larger than that of the OviI line due to SWCX found in the 
north ecliptic pole observation (5.1]';J;o LU, Fujimoto et al. 2007), 

This strong OviI emission line motivated us to investigate the light curve in detail. 
Below we use only the back-illuminated CCD (XISl) data since it has the higher statistics 
for the oxygen line compared to the front-illuminated CCDs as shown in figure 1. The data 
reduction was performed in the same way as Ebisawa et al. (2008) except that the HEAsoft 
analysis package was ver 6.1.1. 

3. Timing Analysis 

To examine the SWCX interpretation, we created the XISl light curve in 0.5-0.65 keV 
where the O VII line dominates and compared it with solar wind data. Similar to Ebisawa et 
al. (2008), we excluded two point sources within the field of view (Suzaku J1844— 0404 and 
G28.6— 1), in order to avoid contamination from these sources. The total area after removing 
two regions is 275 arcmin^. Figure 2 shows the XISl 0.5-0.65 keV light curve in 8192 s bins, 
compared to the WIND proton and ACE O''"'" ion fluxes ^. The position of WIND during the 
observation was -1-200 Re (Earth radius) and —50 i?E in the GSE x-y coordinates (i.e., at the 
pre-bowshock position^), while ACE orbited around the Lagrangian point Li between the Sun 
and the Earth. Because the ACE level 2 (publication-quality) proton data was unavailable 
during a part of the observation, we used the WIND data for the proton flux. The average 
XISl count rate was 1.5 x 10^^ cts s~^. 

We can see an increase of the XISl light curve from 303 to 304 day. We checked the 
XISl light curves in different energy bands (0.65-1 keV and 1-10 keV) but both of them did not 
show such a variability (see fig. 2). Hence, this feature seems to be intrinsic for the OviI line 
emission. A similar enhancement can be seen in the solar wind proton and 0^"^ fluxes. Since 
the O^"*" ion is the source for the SWCX OviI line which is from O^"^, this provides a strong 

^ The WIND and ACE data were taken from ftp://space.niit.edu/pub/plasnia/wind/kp_files/ and 

http: / / www.srl.calteeh.edu/ACE / ASC /level2 / index.html, respectively. 
^ http://edaweb.gsfc.nasa.gov/cgi-bin/gif_walk 
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line of evidence that the O VII emission is truly arising from the SWCX. The heliospheric origin 
can be rejected since the heliosphere is far larger than the Earth's geocorona in size and the 
short time variability of the solar wind will be smeared. Therefore, we conclude that the O VII 
line is likely from the geocoronal SWCX. 

Because the ACE satellite orbits at the Lagrangian point Li and Suzaku is in the low 
Earth orbit, we can expect ~ 1 hr time delay between the XISl and ACE data. In order 
to examine the time delay, we conducted a cross-correlation analysis. This procedure needs 
that both light curves are taken in the equally-spaced time intervals. Since the time bins of 
ACE O''"'' flux and XISl light curve are different, we interpolated the ACE data to match the 
XISl. Considering the ACE O^"*" flux is a 2 hr average, both the XIS and ACE data were 
binned into 8192 s. We then utilized crosscor in the HEAsoft analysis package, to obtain the 
cross correlation. In this software, we can choose several parameters for the calculation of the 
cross correlation function. We chose the default mathematical algorithm (f ast=l, fast Fourier 
transform) and normalization method (normalization=l, no renormalization). 

The calculated cross correlation is shown in figure 3. At the time delay around 0^16384 
s, the correlation coefficient was ~0.72 with the null hypothesis probability of 1 x 10~^ corre- 
sponding to ~ 4cr significance. Hence, the correlations around the time delay of 0^16384 s are 
highly significant. We noticed that the peak is slightly shifted to the positive delay side, which 
means that the ACE data has a time delay against the XIS. This coincides with the fact that 
the ACE satellite orbits at the Li point and hence detects the solar wind before Suzaku. The 
expected time delay roughly depends on the distance between the two satellite (~ 1.5 x 10^ km) 
and the average solar wind proton speed at the observation time (~ 400 km s~^). Then, the 
time delay is estimated as ~3800 s. This is consistent with the observed time delay of 0^16384 
s. The sparse ACE fiux data hindered us to investigate a more accurate determination. Below 
we simply assume the time delay of 8192 s, corresponding to the peak of the cross correlation. 

In figure 4, we plot a relation between the XIS OVII line rate and the ACE O^^ fiux 
considering the time delay of 8192 s. We fitted the data with a linear function (solid line). The 
best-fit function was expressed as, 

Cxisi [cts s"^] = Co-+ [10^ cm-2 s"^] x (7.8 ± 1.3) x 10"^ + (1.0 ± 0.1) x lO'^, (1) 

where Cxisi and Co7+ are the XISl 0.5-0.65 keV count rate and the ACE O''"'' fiux, respectively. 
The error is Icr statistical one. The reduced was 0.32 for 25 degree of freedom and hence 
acceptable. This relation suggests that a linear increase of the XISl count rate according to 
the ACE O^''" fiux, which coincides with the SWCX picture. To know the effect of the assumed 
time delay on the relation, we tested other values (0 and 16384 s). We found that, even in 
these cases, the best-fit linear and constant coefficiencies coincide with those assuming 8192 s 
within 1 a fitting errors. 

We noticed that, even if the solar wind O^"*^ fiux is zero, there remains a certain offset 
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component. This is considered to consist of the instrumental and sky background. To clarify 
the offset component, we estimated the instrumental background using the XIS background 
data base, which was built from night Earth observations (Tawa et al. 2008), in the same way 
as Ebisawa et al. (2008). The XISl 0.5-0.65 keV instrumental background was estimated as 
(3.9 ±0.1) X 10"^ cts s~^ which occupies 40 % of the offset value. The rest must be a stable 
background from the sky, (0.6 ±0.1) x 10"^ cts s~^ corresponding to ~ 4 LU in flux. 

The most plausible interpretation for this sky background is a sum of the heliospheric 
SWCX and the soft X-ray background. The former occurs when the solar wind ions react with 
a neutral target in the heliosphere (e.g., Koutroumpa et al. 2007) but its contribution is unclear, 
since it is observationally difficult to distinguish the heliospheric SWCX from the geocoronal 
SWCX except by time variations. The latter is considered to consist of faint extragalactic 
sources and emission from highly ionized ions in solar neighborhoods and in our Galaxy (e.g., 
McCammon et al. 2002; Masui et al. 2009; Yoshino et al. 2009). Recently, Masui et al. (2009) 
studied the soft X-ray background from the galactic disk (/ = 235°, 6 = 0°) with Suzaku. The 
OVII line flux in this line of sight was 2.93±0.45 LU. The line of sight of our Galactic ridge 
observation (/ = 28.463°, b = —0.204°) is in the direction of the galactic disk. The slightly 
different flux may be due to possible directional dependence of the soft X-ray background 
and/or seasonal or directional differences in the heliospheric SWCX. 

As an additional background, we also considered the fluorescence scattering of solar X- 
rays by the Earth's atmosphere. The fluorescent oxygen line at 0.53 keV is sometimes seen in 
past Suzaku observations (e.g.. Miller et al. 2008). The contribution of the fluorescent emission 
line depends on both the amount of the solar X-rays during the observation and the Earth's 
atmosphere density along the line of sight, which is related to the telescope elevation angle from 
the Earth (ELV). To estimate this influence, we flltered the XISl data with different ELV of 
> 5° (default criteria), > 10°, > 20°, > 30°, and > 50°. We then found no signiflcant change in 
the count rate when the ELV criteria is different. This suggests a negligible contribution from 
the fluorescence scattering. This is consistent with the fact that the solar X-rays were quite 
stable at a low level during the observation (7 x lO^'^ 3 x 10^^ W m^^ in the G0ES12 1.0-8.0 
A data)^. This flux level corresponds to the minimum class solar flare (class A) or less. Hence, 
we concluded that the fluorescent scattering of solar X-rays is negligible in our data. 

4. Discussion 

We analyzed the Suzaku XISl Galactic ridge observation data and discovered a good 
time correlation between the XISl OviI count rate and the ACE O''"'' ion flux. From the line 
center energy and the time correlation with the solar wind, we concluded that this emission 
originates from the geocoronal SWCX. The average OviI line flux was 11 LU, composed of 7 

^ http://www.swpc.noaa.gov/ 
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LU from the geocoronal SWCX and 4 LU from the sky background including the soft X-ray 
background and the hehospheric SWCX. Below, we discuss the expected line intensity of the 
geocoronal SWCX. 

Figure 5 shows the schematic view of the Earth's magnetosphere and line of sight in 
our observation. While the north ecliptic pole observation with Suzaku was in the direction 
of the north pole (Dec= 66°) (Fujimoto et al. 2007), the line of sight in this observation is 
almost vertical to it (Dec= —4°). The line of sight corresponds to the sub-solar side of the 
magnetosheath which is predicted to be a high SWCX flux region by the theoretical model 
(Robertson & Cravens 2003). 

The line of sight of the Galactic ridge also makes a clear difference in the geocentric 
distance of the point whose geomagnetic field becomes open to space for the first time along 
the line of sight (r^p, see fig. 5 for definition). This definition is the same as that in Fujimoto 
et al. (2007). In our case, it is equal to the line-of-sight distance to the magnetopause boundary 
which is the location where the outward magnetic pressure balances the solar wind pressure. 
Here important parameter to estimate the geocoronal SWCX intensity because most 

of the solar wind are blocked by the Earth's magnetic field and hence interact around r^p. 
Since the Earth's exospheric density strongly depends on the distance from the Earth (e.g., 
Ostgaard et al. 2003), r^p should have a large impact on the SWCX intensity. 

Fujimoto et al. (2007) evaluated r^p using the software GEOPACK-2005 and T96 mag- 
netic field model (Tsyganenko et al. 2005) ^. In this model, they took into account solar 
wind parameters, that can influence on the structure of the magnetosphere, using CDAWeb 
(Coordinated Data Analysis Web) ^. We calculated r^p in the same way and found that it is 
almost constant around ~12 Re, corresponding to the distance to the magnetopause. This is 
natural consequence if we consider the line of sight direction (see fig. 5). r^p is larger than 
that in the north ecliptic pole (2~8 Re)- In order to investigate whether the short term vari- 
ability seen in the north ecliptic pole was due to the changing rmp, we checked the short term 
variability less than several hours by changing the bin size of the OVII light curve from 512 s 
to 8096 s in our data. Then we found no sign of such a variability, which is consistent with the 
almost constant r^p. 

We then proceeded to estimate the expected X-ray line intensity from an equation below. 



where /ovii and /07+ are the OviI line intensity and solar wind O^^ flux, rin is the neutral 
geocoronal hydrogen density, and / is the line of sight, /max and /min are the most distant and 
closest distances at which the solar wind ions can interact with the neutral target, a contains 
the atomic cross section and transition probability information. 

^ http: / / niodelweb.gsfc.nasa.gov / magnetos / data-based / modeling. html 
^ littp://cdaweb. gsfc.nasa.gov/cdaweb/sp phys/ 




[photon s ^ cm ^ str ^] 



(2) 
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To calculate /ovii, we made simple assumptions as below. We used a of 6xl0~ cm~ 
from Wegmann et al. (1998), assuming that all transitions are equally probable. The Earth's ex- 
ospheric neutral hydrogen density model by Ostgaard et al. (2003) was used for nii{l). Because 
most of the solar ions will be trapped near the Earth's magtopause in the direction of the 
Galactic ridge, we simply assumed /min = r^p. /max of 20 Re was assumed as a rough estimate, 
because nu{l) at >12 Re is not well known due to the very low density (< 20 cm~'^, see fig. 
10 in Ostgaard et al. 2003). The average ACE O''^ ion flux of ~ 1 x 10~^ cm~^ s~^ was used 
as /o7+- Then, the OVII line flux was estimated as 0.2 LU, which was 35 times smaller than 
the observation. Even if we integrate the equation up to 100 Re, it was 0.3 LU, still 20 times 
smaller. 

In the above calculation, we did not take into account an effect of the bow shock across 
which can increase the solar wind flux. Then, we tried to take into account this effect based on 
a standard theoretical model given by Spreiter et al. (1966). According to their calculation, the 
solar wind mass flux increases 2~3 times across the bow shock in the sub-solar side. Although 
we take into account this effect, there remains a factor of ~10 discrepancy between the observed 
and expected OVII line flux. This uncertainty can be due to an uncertainty of the neutral 
hydrogen model at >10 Re and/or the solar wind distribution. 

The time delay between the X-ray light curve and the solar wind ion flux can provide 
us with information on the propagation of the solar wind. Carter & Sembay (2010) checked 
the ACE and WIND proton curves, and estimated the orientation of the solar wind wavefront 
assuming a planar wavefront. In our case, a part of the ACE proton flux during the time 
variation is unfortunately unavailable. On the other hand, as we described in §3, the time 
delay between the Suzaku light curve and the ACE ion flux of 0~ 16384 sec is consistent with 
the calculation assuming the planar wavefront and constant solar wind velocity. This supports 
that enhancement of the geocoronal SWCX diffuse X-ray background can be predicted and 
excluded if there are simultaneous solar wind proton and ion observations. 

Finally, we point out that the geocoronal SWCX in combination with the solar wind 
data can be a useful tool to evaluate the Earth's exospheric density model and the solar wind 
distribution. Our observation had a line of sight which intersects with the sub-solar side of the 
magnetosheath, which is predicted to have a high SWCX flux. Archival Suzaku data analysis 
will enable us to examine the directional dependence even more. Our observation demonstrated 
that the time delay and intensity relationship between the X-ray O VII line and solar wind O^"*" 
flux must be a key to investigate the hue of sight dependence of the geocoronal SWCX. Future 
high-energy resolution observations with an X-ray microcalorimeter onboard Astro-H will make 
it possible to monitor the conditions of the Earth's exosphere and solar wind more precisely. 

The authors acknowledge the Suzaku XIS instrument and operation team, and the ACE 
SWEPAM/SWICS instrument team and the ACE Science Center. 
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Fig. 1. Background-subtracted XISl (black) and XISO + 2 + 3 (red) spectra of the Galactic ridge region. 
The solid line is the best-fit power-law plus thirteen Gaussian model. The dashed lines are the model 
components. Parameters are described in table 2 of Ebisawa et al. (2008). The bottom panel indicates 
residuals of the data from the model. 
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Fig. 2. XISl light curves in 0.5-0.65 keV and 1-10 keV, solar wind proton flux, and 0^+ flux as a 
function of day of year (DOY) in 2005. The DOY of 301 corresponds to 2005 October 28. X-ray photons 
from calibration sources are excluded for the XISl 1-10 keV light curve. The vertical error bars are at la 
significance. The solar wind proton flux is calculated from WIND SWE data (0.001 day=86.4 sec average), 
while the 0^+ flux is from level 2 ACE SWICS data (2 hr average, corresponding to the minimum time 
bin). Only good data with quality flag were used for the ACE data. 



10 



d 



I 

4— > 
1 

I 



1 r 1 1 1 1 1 1 j 1 1 r 

- 


-1 [ 


f r r [ r 1 1 t 1 1 ■! 1' 














J 1 1 1 1 1 1 i 1 J 1 t 


, t , 


1 

V. 

I 1 » 1 1 1 1 1 1 1 1 1 



-2x10^ 



Time Delay (s) 



10-^ 



2x10^ 
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Fig. 4. Correlation between the XISl O VII line emission and ACE 0^+ flux curves, considering the 8192 
s time delay of the Suzaku data (see text in §3). The vertical error bar is at la significance. The solid 
curve is the best-fit linear function. 
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Fig. 5. Schematic view of tlic magnctospliere and line of siglit of tlie satellite, observed in the GSE x-z 
plane. The Sun angle (Sun-Suzaku-object) during the observation was ^ 67° and the line of sight has a 
+57° azimuth angle, perpendicular to this diagram. 
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